Diffusion-Reorganized Aggregates: Attractors in Diffusion Processes? 
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A process based on particle evaporation, diffusion and re- 
deposition is applied iteratively to a two-dimensional object 
of arbitrary shape. The evolution spontaneously transforms 
the object morphology, converging to branched structures. 
Independently of initial geometry, the structures found af- 
ter long time present fractal geometry with a fractal dimen- 
sion around 1.75. The final morphology, which constantly 
evolves in time, can be considered as the dynamic attractor 
of this evaporation-diffusion-redeposition operator. The en- 
semble of these fractal shapes can be considered to be the 
dynamical equilibrium geometry of a diffusion controlled self- 
transformation process. 

PACS numbers: 05.40.+j, 61.43.Hv, 64.60.Cn 



This letter reports the discovery of a diffusion medi- 
ated process which spontaneously builds a dynamic frac- 
tal equilibrium structure, in contrast with fractal mor- 
phologies linked to far- from- equilibrium processes 
The process is a surface to surface evaporation-diffusion- 
condensation process which conserves the total mass of 
the system. During the time evolution, one observes a 
progressive transformation of the surface through bulk 
diffusion. After a long relaxation period, the system 
reaches a dynamical fractal structure. This structure ap- 
pears as the final equilibrium state towards which any 
initial morphology of M particles will converge after suf- 
ficient evaporation-diffusion-condensation iterations. It 
can then be considered as a general statistical attractor 
for that specific dynamic process. 

The underlying ideas that have suggested this study 
come from our knowledge of the basic mathematical ob- 
jects which govern the exchange of Laplacian driven cur- 
rents across irregular (or fractal) interfaces (as, for exam- 
ple, in the study of irregular or fractal electrodes). This 
problem can be mapped onto the study of the transfer 
of Brownian particles across irregular membranes with 
finite permeability Q]. In this process, Brownian parti- 
cles strike an irregular surface where they are absorbed 
with finite probability. When reflected, the particles un- 
dergo successive random paths, hitting and hitting again 
the static surface, until they are finally absorbed. Halsey 
first indicated that the response of such systems depends 
on the probability that a particle starting on the interface 
comes back to it [||. Generalizing these ideas, it has been 
recently shown that the Laplacian transfer across irreg- 
ular interfaces is controlled by a single linear operator Q 



which maps the static surface onto itself through effective 
"Brownian bridges" |6). In this context, each surface has 
an operator Q, which is symmetric and positive. 

The question arises naturally of the link between the 
Brownian bridges and a real diffusion process. This led us 
to study a dynamical process in which the notion of a Q 
operator, which maps the surface onto itself through av- 
erage diffusion, is transformed into a real operation which 
transforms the surface through discretized diffusion [Q. 
The evolution mechanism then proceeds in the following 
steps (Fig. 1): 

1. A particle is chosen at random on the surface of 
an initial structure and is dissolved, occupying an 
empty site next to its initial position. 

2. The dissolved particle may jump back to its original 
site or start a Brownian random walk on a square 
lattice which represents the outer medium. The 
time step for each jump is r = 1. This motion is 
stopped when the particle hits the structure again 
as it sticks on first hit ||. 

3. The process, which conserves the total mass, is then 
iterated. For practical reasons, the system is en- 
closed in a large square window: whenever a dif- 
fusing particle strikes the window boundary, this 
particle is reinitiated on another site of this exter- 
nal box, following a probability law that simulates 
the random walk in the infinite outer medium. As 
the probability to return to a starting site is equal 
to 1 in d = 2, this does not modify the morphologies 
which are generated by the process. The existence 
of the external box only perturbs the time scale so 
that the number of computer steps is not propor- 
tional to a real time. For this reason, the number 
of steps will be called pseudo-time t. 
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FIG. 1. Schematic of the self- diffusion reorganization pro- 
cess. Step 1: a particle on the boundary of the structure is 
randomly chosen to be dissolved, according to the rule that 
it does not break the connectivity of the structure. Step 2: it 
can go back to its original position or starts a random walk 
on a square lattice. If the walk brings the particle to a surface 
site, the particle sticks to the structure under the condition 
that the new position conserves the simple connectivity of the 
structure. 
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The dissolution-recondensation process rapidly creates 
several branches. If particles of these branches were al- 
lowed to dissolve, it would lead to a progressive splitting 
of the structure into two, then many, disconnected parts. 
In order to avoid this progressive disconnection, one im- 
poses the structure to remain connected throughout the 
process. The only particles allowed to evaporate are then 
the particles pertaining to the surface of the structure 
and that are not "red particles": the red particles are 
defined as the particles which, if eliminated, would dis- 
connect the structure in two or more distincts parts. The 
particles allowed to evaporate are therefore called "blue 
particles". We also impose that an evaporated parti- 
cle cannot stick on a site where it would disconnect the 
outer medium. In fact, one could release this last con- 
straint since it does not modify the long term evolution 
of the structure. It simply avoids to create temporary 
"lakes" in which particles evaporating and redeposited 
would greatly reduce the speed of evolution of the struc- 
ture. Particles on the internal surface of these lakes would 
spend most of the time to go back and forth until the lake 
opens spontaneously. 

The results of the spontaneous evolution of structures 
obtained with different initial morphologies are shown in 
Fig. ^. The left column shows the initial shape and the 
corresponding final morphologies are shown on the right. 
In case (a), the initial morphology is a compact dense 
structure described by a dimensionality equal to 2. In 
case (b), the starting shape is a line with dimensional- 
ity 1. Finally, in case (c), the initial morphology is an 
ordinary Diffusion-Limited Aggregate (DLA), first intro- 
duced by Witten and Sander || . 

In all three cases (a), (b) and (c), the final morpholo- 
gies present branched structures, like DLA but less open. 
The complete statistical study of the morphologies is dif- 
ficult because the process is very time consuming in the 
computer. For instance, the time necessary to reach the 
morphology shown in Fig. |2|(a) is of the order of several 
weeks of CPU time on a Hewlett-Packard CI 60 worksta- 
tion. The similarity in the final shapes, independent of 
initial geometry, suggests that this class of morphologies 
is an attractor for such self-reorganization processes. We 
call these structures "Diffusion- Reorganized Aggregates" 
(or "DRA"). 
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FIG. 2. Three examples of the spontaneous evolution of 
morphologies towards statistical equilibrium shapes under the 
diffusion limited reorganization process. Case (a): the initial 
morphology is dense with dimensionality Df = 2, and con- 
tains 7800 particles. Case (b): initial morphology with di- 
mensionality Df = 1, 983 particles. Case (c): initial DLA 
morphology with dimensionality Df — 1.65, 8000 particles. 
The final morphologies look all the same and keep the same 
statistical characteristics after being formed. The time for 
evolution strongly depends on the initial shape. 

The fact that morphologies are "attracted" by the frac- 
tal final form is also indicated in Fig. ||. The figure gives 
the determination of the fractal dimension of the final 
morphologies starting respectively with initial dimension 
Df = 2,D f = 1, and D f = 1.65 (Fig. | (a), (b) and (c)). 
Although the initial dimension are very different, the final 
plots are within reasonable uncertainty concentrated and 
compatible with Df = 1.740 ±0.02, a value which is sig- 
nificantly different from the fractal dimension of lattice 
animals and compatible with that of the branched 
structures which are the results of a far-from-equilibrium 
process [XX . 
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FIG. 3. Evolution of the fractality. Top: mass to radius re- 
lation for initial morphologies (a), (b) and (c) (a — > o, b — > +, 
c — > *). Bottom: mass to radius relation for the correspond- 
ing stationary states. This plot indicates that, independent 
of the initial dimension, the stationary dimensions are very 
close. 

The dynamics of the restructuration process has sev- 
eral specific properties. To follow the dynamics, one can 
compute two memory functions. First, a density-density 
memory function. For this purpose, the structure is char- 
acterized by a number Si on each site i of the lattice, with 
Si = 1 if the site is occupied by the structure and if 
not. Calling M the mass of the initial structure, which 
is equal to the number of occupied sites, this memory 
function can be defined as: 



C ^ = m £ Si(0) Si{t) 



(1) 



lattice sites i 



This first function does not discriminate between par- 
ticles as Si(0) Si(t) is equal to 1 if a lattice site i is oc- 
cupied both at pseudo-times and t (by any particle) 
and otherwise. The pseudo-time evolution of C(t) is 
shown in Fig. ||. One observes a memory loss which does 
not go to zero. There can be two reasons for incom- 
plete decorrelation. First, it could happen that a part 
of the object is not restructured. Second, it could be 
that the density-density memory does not vanish due to 
an average statistical overlap, even if the structure has 



completely forgotten its initial geometry. To clarify this 
point, one can compute a particle-particle memory func- 
tion which characterizes the speed at which the particles 
of a given structure finally move from their initial posi- 
tions. This second function is 



*(r}(0)-r}(t)) 



(2) 



particles j 



where fj(t) is the position of the particle j at time t. 

The evolution of P(t) is also given in Fig. [|. One ob- 
serves that the particle-particle memory is lost to zero 
very rapidly. When it has reached very small values, one 
can say that the structure has completely lost the mem- 
ory of its microscopic configuration. Apart from these 
two rapid processes, one observes in Fig. f| the very slow 
global evolution. Note that the DRA geometry loses con- 
stantly its memory whereas the DLA growth keeps the 
total memory of its previous structure. 

A more detailed characterization of the rapid regimes is 
displayed in Fig. ||, where the decorrelations are found to 
be exponential, a non-usual time behavior in this context. 



0.8 

CO 

c 
o 

£0.6 

H — 

°0.4 

CD 



0.2 




P(th 




500 1000 1500 2000 2500 
Pseudo-time (Million steps) 

FIG. 4. Time evolution of the memory functions C(t) and 
P(t) during the evolution of the disk towards DRA structure 
(Fig ^a) . These functions are defined by Eq. [l]and Eq. ^ Both 
functions are reset five times during the evolution (positions 
of the peaks). After one million time steps, the correlation 
C(t) is less than 1/2, indicating that more than half of the 
initial points have moved and no more belong to the initial 
structure. 
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FIG. 5. Detailed behavior of the memory functions decays 
C(t) — C(oo) and P(t). In both cases, the final decay is found 
to be exponential, with characteristic times of re ~ 1-1 10 8 
and rp ~ 5.10 6 . This suggests a classical equilibrium situa- 
tion. 

Those results, together with the fractality evolution of 
Fig. ||, indicate that the final DRA morphologies behave 
as stable fixed points. The constantly changing DRA can 
then be viewed as a statistical equilibrium state. In that 
sense, it is analogous to the homogeneous state of an ideal 
gas, where the microscopic interactions between particles 
lead to a statistical macroscopic equilibrium, although it 
is constantly changing its microscopical structure. 

One should also note that the transformation of the Q 
operator, which maps the surface onto itself through av- 
erage diffusion, into a real operation is analogous to the 
dielectric breakdown model of DLA |^2| . In the DBM 
model, the operator, that maps the source at infinity 
onto the growing morphology through average Brown- 
ian paths, is transformed into a real stochastic operation. 
In the model presented herein, the source is the object 
itself. 

It is also interesting to return to the exact concept of 
the self-transport operator attached to a given morphol- 
ogy ||. Extending the notion of diffusive self-transport 
in this context, it might be that DRA shapes could be 
considered as " eigenshapes" of these operators. By this, 
it is suggested here that the self-transport operators at- 
tached to these structures transform those structures into 
themselves, of course in a statistical sense. 

All these results should be confirmed in the future 
through extensive numerical simulations and extension to 
diffusive self-transport in d = 3. Some future extensions 
of this work can be envisaged. First, one should verify 
the universal character of this result. Second, if surface 
energy was included, it would be interesting to study how 
the equilibrium morphology would depend on its value. 
Last, one should study whether the self-transport oper- 
ators Q of these geometries possess some extremal prop- 
erties. 

In summary, it has been shown that diffusive self- 
transport statistically modifies an object morphology in 
an irreversible manner. This irreversible process leads 
to a class of branched structures, called DRA, with frac- 



tal dimension around 1.75. This class of structures plays 
the role of a statistical attractor. Moreover, its attrac- 
tion basin seems to be the whole configuration space in 
the sense that, independent of the initial geometry of the 
system, the dissolution-redeposition process always tends 
to create the same category of geometries. The fractal 
DRA can be seen as a statistical equilibrium for this pro- 
cess, which would maximize a "geometrical entropy" still 
to be defined. 
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